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蝎 虎 天 体 PKS 07354178 的 光 变 特性 分 析 
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摘要 : 在 收集 大 量 数据 的 基础 上 ， 用 时 间 补 偿 离 散 传 里 叶 变 换 、Jurkevich 方法 和 离散 相 
关 分 析 法 分 析 了 PKS 0735+178 的 B 波 段 和 V 波 段 光 变 周期 ， 发 现 该 天 体 具 有 年 的 光 变 周期 ， 
其 中 心 黑洞 质量 的 下 限 为 。 
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BL Lac 天 体 是 活动 星系 核 中 一 个 重要 的 子 类 ,BL Lac 天 体 在 观测 上 表现 出 高 光度 、 高 
偏振 、 快 速 光 变 以 及 非 热 辐 射 等 特征 023, BL Lac 天 体 具 有 长 周期 光 变 和 短 时 标 光 变 B3, 通过 
观测 和 研究 BL Lac 天 体 不 同 的 光 变 时 标 能 获得 天 体 的 中 心 黑 洞 质量 、 辐射 区 域 及 内 部 结 
构 参 数 等 5617。 因 此， 研究 该 天 体 的 光 变 周期 非常 重要 。 

BL Lac 天 体 PKS 0735+178 是 一 颗 红 移 为 0.424 的 泡 变 体 趾 ， 它 具有 平 射电 谱 和 射电 爆 、 
剧烈 的 光 变 、 高 偏振 及 超 光 速 运动 等 性 质押 ， 不 发 出 或 仅 发 出 微弱 和 间断 的 发 射线 09. 众 所 
周知 ， 分 析 PKS 0735+178 光 变 周期 的 方法 很 多 ， 但 是 有 些 方法 误差 极 大 ， 对 观测 数据 的 连 
续 性 有 极 高 的 要 求 ， 并 不 适合 周期 性 研究 ， 而 BL Lac 型 天 体 由 于 观测 特征 表明 无 发 射线 
中 ， 要 通过 光谱 观测 获取 光谱 线 ， 并 用 其 研究 中 心 黑 洞 质量 等 内 部 结构 参数 是 不 可 能 的 2 
3。 本 文 主要 使 用 时 间 补 偿 离散 传 里 时 变换 、 离 散 相 关 分 析 法 和 Jurkevich3 种 方法 对 B 波段 
和 V 波段 的 光 变 周期 进行 研究 ， 并 对 比分 析 这 3 种 方法 ， 其 中 时 间 补 偿 离 散 传 里 时 变换 和 
离散 相关 分 析 法 是 初次 用 于 研究 PKS 0735+178 的 光 变 周期 ， 这 两 种 方法 对 观测 数据 的 连 
续 性 要 求 低 、 结 果 准 确 ， 而 Jurkevich 方法 在 文 [9] 中 用 过 ， 并 且 它 要 求 观 测 数据 的 时 间 序 列 
长 、 连 续 性 强 。 


1 样本 和 光 变 曲线 


本 文 研究 的 数据 从 文 [11-18] 获 取 ， 收 集 从 1970 年 到 2002 年 期 间 PKS 0735+178 光学 B 
波段 和 V 波段 的 观测 数据 。 图 1 和 图 2 分 别 为 B 波段 和 V 波段 的 光 变 曲线 。 从 光 变 曲线 可 
以 看 出 : PKS 0735+178 在 光学 波段 活动 非常 激烈 ， 在 近 30 年 的 观测 中 ，B、V 波段 的 最 大 
t 变化 近 3.5 个 星 等 ， 由 于 受 观测 条 件 的 限制 ， 光 变 曲线 数据 不 连续 ，V 波段 将 近 9 年 没有 数 
- 据 ， 使 周期 性 分 析 受 到 限制 。 
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2 周期 分 析 


2.1 时 间 补 偿 离 散 傅 里 叶 变 换 分 析 PKS 0735+178 的 光 变 周期 
时 间 补 偿 离散 传 里 叶 变换 方法 是 计算 光 变 周期 最 常用 的 方法 之 一 ， 文 [19] 用 该 方法 分 析 
PKS 1510-089 红外 光 变 周期 ,通过 对 、、 作 Gram-Schmidt 正 交 化 ， 得 到 3 个 正 交 向 量 ， 将 
数据 投影 到 3 个 正 交 向 量 上 就 得 到 了 频谱 : 


(1) 
(2) 


(3) 
正 交 化 后 : 


(4) 
(5) 


(6) 
对 于 均匀 采样 的 数据 ， 这 个 过 程 对 应 于 用 三 维 “ 正 弦 + 常 数 ” 模 型 进行 曲线 拟 合 。 如 果 周 
期 大 于 采样 时 间 ， 且 时 间 序 列 足够 长 ， 履 盖 所 有 的 相位 ， 则 有 


a 


功率 谱 强度 为 : 

(8) 
其 中 ,为 离散 傅立叶 变换 的 功率 谱 。 在 不 均匀 采样 的 情况 下 ， 是 加 权 的 时 间 补 偿 离散 传 立 叶 
变换 : 


(9) 
在 许多 蝎 虎 天 体 的 观测 中 ， 观 测 数据 的 精度 各 不 相同 。 考 虑 到 这 个 问题 ， 引 进 权重 方程 重 
新 定义 内 积 : 


(10) 
频率 处 的 强度 由 下 式 给 出 : 


(11) 
由 线性 回归 理论 可 知 , 其 中 。 利 用 这 一 性 质 ， 引 进 标准 化 因子 : 统计 量 , 称 这 一 量 为 谱 相 关 
系数 ， 对 于 所 有 频率 ，。 


DCDFT(B) 
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应 用 上 述 方法 获得 PKS 0735+178 的 B 波 段 和 V 波段 周期 图 ， 由 图 3、 图 4 可 知 ，B 波 
段 存在 0.62 年 、1.01 年 、1.67 年 、4.72 年 的 光 变 周期 信息 ，V 波段 存在 0.58 年 、1.07 年 、 
1.67 年 4.72 年 的 光 变 周期 信息 ， 按 文 [19] 中 的 判 据 ， B 波段 存在 1.01 EI V 波段 存在 1.07 
年 光 变 周期 信息 的 可 靠 性 最 大 , 而 B 和 V 波段 存在 4,72 年 次 之 .但 考虑 到 平均 时 标 为 年 的 周 
期 很 可 能 是 由 于 地 球 绕 太 阳 所 致 , 则 光 变 周期 B 波段 为 4.72 年 ，V 波段 为 4.72 年 . 

2.2 Jurkevich 分 析 周 期 

Jurkevich 基于 对 天 文 测 量 中 的 非 均 匀 测 量 问 题 提 出 的 一 种 统计 法 外， 如 果 观 测 样 本 数 
据 为 N 个 ， 为 单 次 测量 值 ， 为 所 有 测量 值 的 平均 值 ， 为 测量 数据 的 方差 ， 为 测量 数据 样本 
的 标准 偏差 . 则 有 


” (12) 


” (13) 


. (14) 
如 果 样 本 数据 划分 为 m 组 ， 对 应 第 1 组 的 统计 参数 为 


' (15) 
' (16) 


. (17) 
对 应 m 组 的 总 方差 为 : 


(18) 


0 , n 1 | 1 1 1 1 1 0 n i i 1 1 1 1 1 1 
0 200 400 600 800 1000 1200 1400 1600 1800 2000 0 200 400 600 800 1000 1200 1400 1600 1800 2000 
P(days) P(days) 


应 用 上 述 Jurkevich 方法 获得 PKS 0735-178 的 B 和 V 波段 -图 ， 由 图 5、 图 6 可 以 看 
出 ， 对 应 于 的 最 小 值 ，B 波段 存在 1.02 年 、1.20 年 、1.66 年 、4.77 年 的 光 变 周期 ，V 波段 
周期 存在 0.94 年 、1.71 年 、2.82 年 、3.97 年 的 光 变 周期 。 根 据 光 变 周 期 存在 的 判 据 中 ，B 
波段 存在 1.02 年 和 V 波段 存在 0.94 年 光 变 周期 信息 的 可 靠 性 最 大 ,而 B 波段 存在 4.77 和 V 
波段 存在 3.97 年 次 之 .但 考虑 到 平均 时 标 为 年 的 周期 很 可 能 是 由 于 地 球 绕 太 阳 所 致 , 则 B 波 
段 存在 4.77 年 周期 ，V 波段 存在 3.97 年 周期 。 

2.3 离散 相关 分 析 法 分 析 
离散 相关 分 析 法 是 用 来 分 析 两 组 离散 数据 相关 性 的 方法 之 一 外 3, 文 [24] 用 此 方法 分 析 BL 
Lac 天 体 PKS 0537-441 的 光 变 特性 史 ， 该 方法 最 大 的 特点 是 不 需要 对 数据 做 任何 处 理 ， 就 
能 判断 出 两 组 数据 的 相关 性 。 该 方法 的 具体 的 步骤 如 下 : 
首先 计算 两 组 数据 的 离散 相关 函数 值 。 如 数组 和 ， 则 离散 相关 函数 值 为 : 
. (19) 
和 分 别 为 两 组 数据 的 平均 值 ， 分 别 为 对 应 的 标准 偏差 。 

其 次 计算 DCF(t) 值 。 通 过 时 间 延 迟 把 两 组 数 联系 起 来 ， 假 如 时 间 延 迟 为 r+， 在 区 间 中 有 M 


chinaXiv:201805.00173v1 


ChinaXiv& 


^r. 则 DCEO 7S 


. (200 
再 次 离散 相关 函数 的 误差 为 
. (21) 


对 于 所 得 到 的 离散 相关 图 ， 如 果 峰 值 在 零 的 右面 ， 表 明 数 组 早 于 数组 的 变化 。 反 之 ， 数 组 
就 会 迟 于 数组 。 


-3 1 1 
-1500 -1000 -500 


0 500 1000 1500 delay(day) 
delay(day) 


通过 使 用 离散 相关 函数 对 PKS 07354178 的 B 波段 图 像 进行 分 析 ， 所 得 结果 如 图 
7、8， 从 图 中 得 到 B 波段 周期 存在 1.21 年、 2.27 年 、3.62 年 、3.81 年 的 光 变 周期 信息 ，V 
波段 存在 周期 1.21 年 、2.22 年 、3.15 年 、4.00 年 的 光 变 周期 信息 。 根 据 离散 相关 示 数 分 析 
周期 的 方法 2 B V. 波段 存在 1.21 年 光 变 周期 信息 的 可 靠 性 最 大 ,而 B 波段 存在 4.77 和 
V 波段 存在 3.97 年 次 之 .但 考虑 到 平均 时 标 为 年 的 周期 很 可 能 是 由 于 地 球 绕 太 阳 所 致 , 则 
PKS 0735+178B 波段 的 周期 为 3.81 年 ，V 波段 的 周期 为 4.00 F- 


3 讨论 与 结论 


作 期 刊 


本 文 主要 收集 了 PKS 0735+178 的 光 变 数据 ， 运 用 时 间 补 偿 离 散 侍 里 叶 法 变换 、 


urkevich 方法 和 离散 相关 分 析 法 3 种 方法 分 析 PKS 0735+178 的 光 变 周期 ， 时 间 补 偿 离散 
傅 里 时 法 变换 法 分 析 也 波段 得 到 的 周期 大 致 为 ,V 波段 周期 为 ，Jurkevich 方法 分 析 得 到 B 波 
段 的 周期 大 致 为 1740 天 ，V 波段 的 周期 大 致 为 1450 K; DCF 方法 分 析 的 得 到 B 波段 的 周 
期 大 致 为 1390 天，V 波段 的 周期 大 致 为 1460 天 .由 此 可 得 出 ，BL Lac 天 体 PKS 0735+178 
的 光 变 周期 为 。 这 与 文 [25] 和 文 [11] 用 其 他 方法 研究 结果 基本 一 致 。 根 据 周期 分 析 获 得 的 
PKS 0735+17 的 周期 ， 我 们 利用 薄 吸 积 盘 理论 B 来 分 析 该 蝎 虎 天 体 PKS 0735+178 的 中 心 黑 
洞 质 量 和 注 吸 积 盘 的 热 不 稳定 性 发 生 的 区 域 。 一 般 情况 下 ， 热 有 限 循 环 周期 性 取决 于 烙 滞 
度 参 数 ， 中 心 黑洞 质量 为 和 广义 应 力 张 量 参数 ， 因 此 爆发 时 间 为 
22 

其 中 ,为 粘 清 系 数 ;， 以 倍 太 阳 质 量 为 单位 ;为 PKS 0735+178 的 中 心 黑洞 质量 。 

由 于 目前 吸 积 盘 粘 滞 度 的 起 源 和 特性 都 还 不 清楚 ， 用 磁 流 体力 学 对 其 讨论 是 一 种 常用 
的 方法 ， 文 [26] 提 出 如 果 磁 场 的 逃逸 率 比 较 低 ， 可 以 认为 ， 使 用 这 个 参数 值 ， 长 周期 爆发 
时 间 大 约 为 ， 即 


(23) 
对 于 PKS 0735+178， 如 果 ，， 分 析 获 得 的 周期 是 ， 可 以 得 出 其 中 心 黑洞 质量 为 ， 这 个 中 心 
黑洞 质量 对 于 PKS 0735-178 来 说 太 小 ， 由 于 利用 薄 吸 积 盘 理 论 分 析 方 法 考虑 黑洞 的 自 旋 ， 
所 以 中 心 黑洞 质量 偏 小 ， 但 蝎 虎 天 体 的 观测 特性 表明 ， 这 类 天 体 在 光谱 观测 中 无 发 射线 
09， 作 为 深入 研究 该 类 天 体 给 出 了 黑洞 质量 的 下 限 也 是 有 意义 的 ， 在 今后 的 CCD 测 光 观测 
中 寻找 短 时 标 光 变 ， 从 而 获得 中 心 黑洞 质量 ， 与 之 做 进一步 比较 。 


从 图 上 也 可 看 出 不 同 的 方法 周期 性 的 明显 程度 也 不 一 样 。 时 间 补 偿 离 散 傅 里 叶 变 换 法 可 以 
运用 到 天 体 数 据 处 理 上 ， 可 以 有 效 地 处 理 不 均匀 的 天 文 数据 ， 克 服 传 统 的 全 里 叶 变 换 带 来 
的 伪 周 期 ， 时 间 补 偿 离 散 传 里 叶 变换 法 通过 施 密 特 正 交 化 ， 有 效 地 解决 由 于 、 和 3 个 向 量 
不 正 交 而 带 来 伪 信 号 的 问题 四 。 用 离散 相关 性 函数 对 PKS 0735+178 天 体 的 B 波 段 和 V 波 
段 相关 性 分 析 ， 从 图 中 可 以 看 出 它们 之 间 存 在 很 强 的 相关 性 ， 表 明 它 们 之 间 的 辐射 起 源 是 
相似 Ie, 
另外 ， 在 光 变 周期 信息 分 析 中 ， 注 意 到 年 这 个 结果 可 能 和 一 年 的 观测 周期 有 关 ， 这 个 周期 
可 能 是 由 于 有 规律 的 观测 时 间 间 隔 导 致 的 中， 平均 时 标 为 年 的 周期 很 可 能 是 由 于 地 球 绕 太 
所 致 中 。 当 然 ， 也 可 能 是 天 体 PKS 0735+178 的 真实 光 变 周期 ， 其 他 光 变 周期 信息 是 这 个 
周期 的 N fii CN=1，2，3...) 四。 这 个 周期 的 研究 四， 对 今后 的 观测 有 指导 意义 。 将 通过 
实测 和 更 多 的 观测 证 明天 体 PKS 0735+178 的 光 变 周期 。 
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Abstract:In this paper, we used the time compensation discrete Fourier transform method and 
discrete correlation analysis (DCDFT), Jurkevich ,DCF model to analyses B band and V band 
periodicity of the PKS 0735 + 178, based on a large amount of data we have collected. We found 
that the variability period is about years and the lower limits of the central black hole mass is 
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